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On the origin of fluorine in the Milky Way 
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ABSTRACT 

The main astrophysical factories of fluorine ( 19 F) are thought to be Type II super- 
novae, Wolf-Rayet stars, and the asymptotic giant branch (AGB) of intermediate 
mass stars. We present a model for the chemical evolution of fluorine in the Milky 
Way using a semi-analytic multi-zone chemical evolution model. For the first time, 
we demonstrate quantitatively the impact of fluorine nucleosynthesis in Wolf-Rayet 
and AGB stars. The inclusion of these latter two fluorine production sites provides a 
possible solution to the long-standing discrepancy between model predictions and the 
fluorine abundances observed in Milky Way giants. Finally, fluorine is discussed as a 
possible probe of the role of supernovae and intermediate mass stars in the chemical 
evolution history of the globular cluster ui Centauri. 
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1 INTRODUCTION 

The three primary astrophysical factories for fluorine ( 19 F) 
production have long been thought to be Type II Super- 
novae (SNe II), Wolf-Rayet (W R) stars, and asymptotic 
giant branch (AGB) sta r s (e.g. IWooslev fc Weaver! Il995t 



iMevnet fc Arnouldl 1200(1 : iForestini et alJ Il992l : Mowlavi, 
Jorissen & Arnould 1998, respectively). Previous attempts 
to model the Galactic production and evolution of 19 F have 
been restricted to explore the role of SNe II alone (e.g. 
Timmes, Woosley & Weaver 1995; Alibes, Labay & Canal 
2001). 

The above problem has now been ameliorated by the 
release of the first detailed yield predictions for fluorine pro- 
duction from WR and AGB stars. We are now in a position 
to incorporate these yields into a Galactic chemical evolu- 
tion framework, in order to assess the respective contribu- 
tions of the three putative fluorine production sites. To do 
so, we will make use of GEtool, a semi-analytical multi-zone 
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Galactic chemical evolution package which has been cali- 
brated with extant observational data for the Milky Way 
jFenner fc Gibsor]l2005t iGibson et al1l2003l) . 



Specifically, in what follows, we compare the model 
fluorine distribution in the Milky Way with the abun- 
dances observed by Jorissen, Smith & Lambert (1992) in 
near-solar metallicity giants. Further, our model predic- 
tions are contrasted with new fluorine determinations for 
giants in the Large Mag ellanic Cloud (LMC) and ui Cen- 
tauri llCunha et al.l l2003). In additi on, ne w results for more 
u Centauri giants from lSmith et alJ i2004T) are included. The 
latter two systems are likely to have had very different star 
formation and chemical evolution histories from those of the 
Milky Way, but despite these obvious differences, a compar- 
ison against these new data can be valuable. In Section 2, 
we provide a cursory overview of the three traditional 19 F 
nucleosynthesis sites; the chemical evolution code in which 
the nucleosynthesis products from these factories have been 
implemented is described in Section 3. Our results are then 
presented and summarised in Sections 4 and 5, respectively. 
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2 NUCLEOSYNTHESIS OF 19 F 

2.1 Type II Supernovae 

The massive star progenitors to SNe II produce fluorine 
primarily as the result of spallat ion of 20 Ne by n and r 
neutrinos near the c ollapsed core (IWooslev fc Haxtonll98St 
IWooslev et aDll99(^ . A fraction of the 19 F thus created is 
destroyed by the subsequent shock but most is returned 
to the ambient Interstellar Medium (ISM). The fluorine 
yields by neutrino spallation are ver y sensitive to the as - 
sumed spectra of fi and r neutrinos iWooslev et al.ll2002j) . 
which could be nonthermal and deficient on their high- 
energy tails, lowering the equiv alent temperature of the neu- 
trinos in the supernova model llMvra fc Burrowslll99Cft . An 
additional source of 19 F derives from pre-explosive CNO 
burning in helium shell. However, fluorine production by 
neutrino spallation is largely dominant, as evide nt by com- 
paring the models in IWooslev fc Weaver! l) 1995(1 , the only 
ones to-date including neutrino process, and recent mod- 
els which do not include neutrino nucleo synthesis of fluorine 
jLimongi fc Chieffil20otl . Most recently. iHeeer et al.l(l2004f) 
suggest that the relevant neutrino cross sections need to be 
revised downwards; if confirmed, the associated SNe II 19 F 
yield would decrease by ~ 50%. In light of the preliminary 
nature of the Heger et al. claim, we retain the conservative 
choice offered by the lWooslev fc Weaverl lll995l) compilation. 

2.2 Asymptotic giant branch stars 

The nucleosynthesis pathways for fluorine production within 
AGB stars involve both helium bu rning and combined 
hydrogen- helium bu r ning phases (e.g . iFor estini et all Il992t 
IJorissen et al] Il992t iMowlavi et al] 1^98^ and are com- 
panions for the nucleosynthesis by s low neutron accre- 
tion (s-process) iMowlavi et al] [l998). Provided a suit- 
able source of protons is available, fluorine can be synthe- 
sised via 14 N(q,7) 18 F (j3 + ) ls O(p,a) 15 N(a, 7 ) 19 F. Primary 
sources of uncertainty in predicting fluorine nucleosynthe- 
sis in AGB stars relate to the adopted reaction rates, es- 
pecially 14 C(a,7) ls O and 19 F(a, p) 22 Ne, and the treatment 
of the nucleosynthesis occurring during the convective ther- 
mal pulses. Nucleosynthesis during the interpulse periods 
can also be important if protons from the envelope are par- 
tially mixed in t he top layers of the He intershell (partial 
mixing zone), as lLugaro et, al.l l|2004h have recently demon- 
strated. Nucleosynthesis in this zone may result in a signifi- 
cant increase in the predicted 19 F yields. The magnitude of 
these systematic uncertainties for stellar models with mass 
~ 3 M and metallicities Z = 0.004-0.02 are ~ 50%, while 
for stellar models with mass M = 5 Mq and metallicity Z 
= 0.02 the uncertainty is a factor of ~ 5, due to the uncer- 
tain 19 F(a,p) 22 Ne reaction rate. Characterising the mass- 
and metallicity-dependence of the partial mixing zone- 19 F 
relationship needs to be completed before we can assess 
its behaviour self-consistently within our chemical evolution 
model of the Milky Way. For the present study, we have 
adopted the yields presented in Appendix, based upon the 
Karakas & Lattanzio (2003, and references therein) mod- 
els, which themselves do not include 19 F nucleosynthesis via 
partial mixing. This choice is a conservative one, and thus 
should be considered as a lower limit to the production of 
19 F from AGB stars. 



For stars more massive than w 4 Mq, the convective 
envelope is so deep that it penetrates into the top of the 
hydrogen-burning shell so that nucleosynthesis actually oc- 
curs in the envelope of the star. Such "hot-bottom-burning" 
acts to destroy 19 F, and should be treated self-consistently 
within the AGB models considered. 

2.3 Wolf-Rayet stars 

Fluorine production in WR stars is tied to its nucleosyn- 
thesis during the helium-burning phase. At the end of this 
phase though, significant fluorine destruction occurs via 
19 F(a,p) 22 Ne. Any earlier synthesised 19 F must be removed 
from the stellar interior in order to avoid destruction. For 
massive stars to be significant contributors to net fluorine 
production, they must experience mass loss on a timescale 
that allows the removal of F before its destruction. This 
requirement is met by WR stars. 

Recently, iMevnet fc Arnouldl J2000l) studied the role 
that such stars can play in the chemical evolution of fluorine 
by adopting updated reaction rates coupled with extreme 
mass-loss rates in not-rotating stellar models. They pointed 
out that WR mass-loss is strongly metallicity-dependent, 
and that the number of WR stars at low metallicities is very 
small. Their WR yields reflect such metallicity-dependence, 
with minimal fluorine returned to the ISM at low metallic- 
ities, but significant 19 F returned at solar and super-solar 
metallicities. The WR yields are sensitive to the adopted 
reaction and mass-loss rates, while rotating models could 
favour an early entrance into the WR phase for a given mass, 
decrease the minimum initial mass for a star to go through a 
WR phase at a given metallicity, and open more nucleosyn- 
thetic channels because of the mixing ind uced by rotation. 
Therefore, after IMevnet fc Arnouldl i2000f) . we consider the 
aforementioned WR yields as lower limits. 



3 THE MODEL 

In this study we employ GEtool, our semi-analytical multi- 
zone chemical evoluti on package to model a sample M ilky 
Way- like disk galaxy jFenner fc Gibsonll2003t iGibson et al.l 
2003). A dual-infall framework is constructed in which the 
first infall episode corresponds to the formation of the halo, 
and the second to the inside-out formation of the disk. 

A Kroupa, Tout & Gilmore (1993) initial mass function 
(IMF) has been assumed, with lower and upper mass limits 
of 0.08 Mq and 120 Mq, respectively. Stellar yields are one 
of the most important features in galactic chemical evolution 
models, yet questions remain concerning the precise compo- 
sition of stellar ejecta, due to the uncertain role played by 
processes including mass loss, rotation, fall-back, and the 
location of the mass cut, which separates the remnant from 
the ejected material in S Ne. The SNe II yields are from 
IWooslev fc Weaverl il995l) : the yields for stars more massive 
than 60 Mq are assumed to be mass-independent. Such as- 
sumption is made to avoid extreme extrapolation from the 
most massive star in the Woosley & Weaver models (40 Mq) 
to the upper end of the IMF (120 Mq), and has negligible 
effect on the results, given the shape of the adopted IMF. 

We have h a lved the iron yield s shown in 
IWooslev fc Weaverl dl995ll , as suggested by iTimmes et al] 
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Figure 1. Fluorine yields from a) SNe II (Woosley & Weaver 
1995), b) WR (Meynet & Arnould 2000), and c) AGB stars (Ap- 
pendix). 



Jl995l) . The Type la (SNe la) yields of llwamoto et alJ jl999l) 
were als o employed. We adopt ed the metallicity-dependent 
yields of lRenzini fc Voli ( 198l|) for single stars in the mass 
range 1 — 8 M©. For the purposes of this work, which 
focuses on fluorine, the choice of the Renzini & Voli yield 
set does not affect the results. Metalllicitv-dependen t stellar 
lifetimes have been taken from ISchaller et al.l (Il992l) . 

We have constructed three Milky Way (MW) model 
variants that differ only in their respective treatments of 
19 F production: 1) MWa assumes that SNe II are the only 
source of 19 F; 2) MWb includes yields from both SNe II and 
WR stars; 3) MWc includes all three sources of fluorine - 
SNe II, WR, and AGB stars. 

We end by noting that within our adopted dual-infall 
framework for the chemical evolution of the Milky Way, our 
model is constrained by an array of observational boundary 
conditions, including the present-day star and gas distri- 
butions (both in density and metallicity) , abundance ratio 
patterns, age- metal licity relation, and G-dwarf distribution 
jGibson et aTl feocm. While the modification of any individ- 
ual ingredient within model framework will have an impact, 
to some degree, upon the predicted chemical evolution, this 
can only eventuate at the expense of one or more of the afore- 
mentioned boundary conditions that we require our model 
to adhere to. Within our framework, yield uncertainties will 
dominate the systematic uncertainties for the predicted evo- 
lution of 19 F. 



3.1 Fluorine yields 

We now summarise the 19 F yields employed in our three 
"Milky Way" models. 

1) SNe II 19 F yields are taken from lWooslev fc Weaver! 



(1995) and assumed to be mass-independent for stellar 
masses in excess of 60 Mq. 

2) WR 19 F yields are taken from iMevnet fc Arnouldl 
(2000) for stellar masses in the range 25 — 120 Mq: each 
star within this range is assumed to evolve through the WR 
stage. Such simplifying assumption could overestimate the 
WR contribution to fluorine, even though the adopted WR 
yields are themselves lower limits (Section 2.3). The WR 
fluorine contribution has been added to the corresponding 
SNe II contribution (which comes from a different stage of 
the stellar evolution). 

3) AGB 19 F and oxygen yields in the 1 - 6.5 M 
mass range have been derived from stellar models con- 
structed with the Mount Stromlo Stellar St ructure Code 
jFrost fc Lattanzidll996l : iKarakas et ai1l2002h , and are pre- 
sented in Appendix. The post-processing nucleosynthesis 
models with 74 species and time-depen dent diffusive con- 
vect ive mixing are described in detail in lFrost et al.l <ll998l) 
and lKarakas fc Lattanziol i2003l) . 

To ensure internal consistency, we have also employed 
the AGB oxygen yields in lieu of those of Renzini & Voli 
(1981), within this mass range. 

The above fluorine yields are shown in Figure 1. In Fig- 
ure 2, the yields are expressed as [F/O] and ([F/0])/a/f, 
the latter corresponding to the mean [F/O] yields for SNe II 
and AGB stars, weighted by the IMF over the SNe II and 
AGB mass range, respectively. We have not shown a compa- 
rable entry for the WR stars as a self-consistent treatment 
of the oxygen production was not included in Meynet & 
Arnould (2000). Here, oxygen has been used as the normali- 
sation to make easier the comparison with the observations, 
especially in uj Centauri, though oxygen can be synthesised 
in various stellar sites, and its yields can be affected by dif- 
ferent reaction rates and modeling of helium cores, semi- 
convection, convective boundary laye rs, and mass-loss (e.g. 
IWooslev et ai]l2002l : iDrav et alJEioa) . 



4 RESULTS 

In Figure 3, the evolution of [F/O], A(O), the gas infall 
rate 6"; n / a ;i, the star formation rate (SFR), the SNe II rate 
and the gas-phase global metallicity Z of the three models 
at the solar neighbour hood are summarised. The empirical 
SFR history derived bv lBertelli fc Nasi <l200ll) is shown as a 
thick solid line in Figure 3d, while the sha ded reg i on co rre- 
sponds to the range of values suggested bv lRanal(ll99ll) . A 
conservative range of estimated SNe II rates is also shown 
in Figure 3e (Cappellaro, Evans & Turatto 1999). 2 . Figure 4 
then shows the the evolution of [F/O] versus A(O) {panel a), 



1 Hereafter, [X/Y]= Zogi O (X/Y)-Zogi O (X/Y) and A(X)= 12 + 
iogio(nx/nH). An accurate determination of photospheric solar 
abundances requires detailed modeling of the solar granulation 
and accounting for departures from local thermodynamical equi- 
librium (e.g. Allende Prieto, Lambert & Asplund 2001). We adopt 
the solar fluorine abundance suggested by Cunha et al. (2003), 
and the solar iron and oxygen abundances from Holweger (2001). 

2 The range of values shown in Figure 3e is derived from the 
sample of SOa - Sb galaxies in Cappellaro et al. (1999) - 0.42 ± 
0.19 SNu, where 1 SNu= 1 SN(100 yr)- 1 (10 10 Lg)" 1 ), assuming 

2 X 10"" 



T B 



10 Lg and a galactic radial extent of 15 kpc. 
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Figure 2. [F/O] and ([F/0]) IMF for SNe II and AGB yields 
(upper and lower panels, respectively). Here, A( 19 F)q = 4.55 
(see discussion in Cunha et al. 2003) and A(0)q = 8.736 (e.g. 
Holweger 2001). The shaded regions in Figures 2a and 2b show 
the observed [F/O] in u) Cen giants (Cunha et al. 2003). The 
([F/0])jMf are weighted by the IMF over the SNe II (11-40M Q ) 
and AGB (1 — 6.5 Mq) mass range, respectively. In Figure 2d, 
both ([F/0]agb)imf and ([F/0]sjv e ii)imf are shown (closed 
boxes and open triangles, respectively). 



Figure 3. Predicted evolution in the solar neighbourhood of a) 
the gas infall rate o"; n / a ;ij b) A(O), c) [F/O], d) star formation 
rate (SFR), e) SNe II rate, and f) metallicity Z (MWa, solid line; 
MWb, dotted; MWc, short-dashed). The SFR history at the solar 
neighbourhood obtained by Bertelli & Nasi (2001) is also shown 
as a thick solid line in panel 'd', while the shaded region shows 
the range of values suggested by Rana (1991). A range of values 
corresponding to the estimated SNe II rate is shown in panel 'e' 
(Cappellaro et al. 1999). 



and the evolution of [F/O] versus [O/Fe] (panel b), compared 
against the IMF- weighted SNe II yields (recall Figure 2). 

The MWa model provides a satisfactory reproduction of 
the estimated star formation history and SNe II rate in the 
solar neighbourhood (Figures 3d and 3e) . This model, whose 
only fluorine source is SNe II, underproduces fluorine with 
respect to the abu ndances measur e d in K and M Milky Way 
giants observed bylJorissen et all <ll992T) and reanalysed by 
ICunha et alJ i2003t) (Fig. 4a). Fig. 4a does not show the s- 
pr ocess enriched AGB stars of spectral types MS, S, or C 
in ljorissen et all lll992l) . where freshly synthesised fluorine 
could be mixed to the stellar surface. Such inclusion of self- 
polluted 19 F-rich stars could obscure any metallicity trend. 
The results of the MWa and MWb models show that the 
additional contribution from WR stars increases [F/O] by 
up to factor of 2 by the present-day, but it is negligible in 
excess of ~ 9 Gyr ago (Figure 3c). 

The addition of both WR and AGB sources within the 
MWc model leads to a present-day [F/O] that is ~ 0.4 dex 
greater than in the MWa case. Further, and perhaps more 
important, AGB stars are now shown to deliver significant 
amounts of fluorine to the ISM during the early epochs of 
the Milky Way's evolution. Such a result is entirely con- 
sistent (and expected) given the metallicity-dependence of 

Given these assumptions, the estimated SNe II rate at the solar 
neighbourhood is necessarily uncertain. 



the AGB yields; said yields possess [F/O] ratios which are 
greater at lower metallicities (recall Figure 2). We can con- 
clude that it is only the addition of both the WR and AGB 
contributions which allow for a significant improvement in 
the comparison between galactic models incorporating fluo- 
rine evolution and the observational data. 



5 DISCUSSION 

We have studied the Galactic chemical evolution of fluo- 
rine, for the first time using new grids of stellar models 
which provide self-consistent predictions of fluorine nucle- 
osynthesis for stars in both the WR and AGB phases of 
stellar evolution. We have shown that the WR contribu- 
tion is significant at solar and super-solar metallicities be- 
cause of the adopted metallicity-dependent mass-loss pre- 
scription employed in the stellar models. In contrast, the 
contribution of AGB stars to fluorine production peaks dur- 
ing the early epochs of the Galaxy's evolution (again due to 
the metallicity-dependent behaviour of the AGB models). 
In combination, the addition of the WR and AGB yields 
leads to a significant improvement in the galactic chemical 
evolution models when compared against observations. 

The comparison between our MW model s and the flu- 
orine abundances in LMC and u> Cen giants iCunha et alJ 
2003) is not straightforward, as the latter two have star for- 
mation (and therefore chemical evolution) histories different 
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Figure 4. (a): [F/O] as a function of A(O) for the MW mod- 
els (MWa, solid line; MWb, dotted; MWc, short-dashed). Also 
shown are the values observed in Milky Way, LMC, and lu Cen 
giants (crosses, boxes and hexagons, respectively), (b): [F/O] as 
a function of [O/Fe], compared with the IMF- weighted [O/Fc] 
yields for SNe II (open triangles). Within the open triangles, '0' 
corresponds to Z=0; '1', to Z=1.9xl0" 6 ; '2', to Z=1.9xl0~ 4 ; '3', 
to Z=1.9xl0 -3 ; '4', to Z=1.9xl0 -2 . The upper panels represent 
enlargements of the framed regions delineated in the correspond- 
ing bottom panels. 



from that of the MW. However, it is interesting to speculate 
on the possible origin of fluorine in u) Cen, given th e unique 
nature of this "globular cluster" (e.g. ISmithll2003|) . Specifi- 
cally, u Cen is the most massive Galactic cluster, and unlike 
most globulars, possesses a significant spread in metallicity 
(~ 1.5 dex) amongst its stellar population. It has been sug- 
gested that u) Cen is ac tually the remnant core of a tidally- 
disrupted dwarf galaxy fBek ki fe Freemanl2003l) . Such a sce- 
nario could naturally drive radial gas inflows to the dwarf 
nucleus, potentially triggering starbursts. 

Interestingly, SNe II ejecta are characterised by low 
([F/0])imf (Figure 2c) and high <[0/Fe]) /A /F (Figure 4b), 
whereas AGB ejecta have higher ([F/OJ);mf (Figure 2d). 
The observed w Cen giants have primarily low [F/O] (Fig- 
ures 2a - 2b, 4a - 4b) and high [O/Fe] (Figure 4b) values, 
consistent with a picture in which their interiors have been 
polluted by the ejecta of an earlier generation of SNe II, but 
not from a comparable generation of AGB. Given that the 
observed oxygen abundance in such u> Cen giants is simi- 
lar to that seen in comparable LMC and MW giants (Fig- 
ure 4a), this would suggest that the chemical enrichment of 
u> Cen proceeded on a short timescale (to avoid pollution 
from the lower mass progenitors to the AGB stars) and in 
an inhomogeneous manner (given the significant scatter in 
observed fluorine a bundances), as discussed previously by 
ICunha et alJ J2Q03). Should the (downward) revised neu- 
trino cross sections alluded to in Section 2.1 be confirmed 



(Heger et al. 2004), the concurrent factor of ~ 2 reduction 
in SNe II 19 F production would improve the agreement of 
the model with the observed F/O ratio in ui Cen giants. 
This would consequently strengthen our conclusions which 
already support a picture whereby these giants have been 
polluted by earlier generations of SNe II ejecta. 
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APPENDIX A: ASYMPTOTIC GIANT 
BRANCH 19 F AND 16 17 lg O YIELDS 

The yields employed here have been derived via the follow- 
ing: 

Yw F> o(Z) = Y netWp ^ o + 

X I19F , o(^) X t m * ~ m * rem(Z)]. (Al) 

Here, Yi9 F: o(Z) is the overall yield, Xi lgp a (Z) is the ini- 
tial mass fraction of the element within a star of mass m* 
and metallicity Z, [m* — r em(Z)] is the total mass ejected 
during the stellar lifetime, and Y ne t WF Q is the net yield. We 
calculate Xi a and Y ne t as, respectively: 

X io = X il6o + Xi 17o + X ilSo ; (A2) 
Y ne t = Y ne t 16o + Y ne t 17o + Y ne t 1So . (A3) 

The yields are the result of full evolutionary calcula- 
tions using the Mount S t romlo Stellar Structure Code 
(e.g. iKarakas fe Lattanziol 12003ft . We use the standard 
Rei mers mass-loss formula o n the first giant branch and 
the lVassiliadis fe Wood! lll993T) for mula during the AGB evo - 
lution. Opacities are from OPAL ^ csias fc Rogers! Hill . 
The models with Z = 0.02 and 0.0001 used scaled solar abun- 
dances, whereas those for Z = 0.004 and 0.008 are appro- 
priate to the Small and Large Magellanic Cloud s, respec- 
tively, and are taken from lRussell fc Dopital <ll992l) . Numer- 
ical problems du ring the third dredge-up are handled in the 
way described in lFrost fc Lattanziol lll996l) . A mixing length 
of 1.75 pressure scale-heights has been used. A main uncer- 
tainty in the predicted yields for fluorine is the occurrence 
and dimension of the partial mixing zone. Note that this par- 
tial mixing zone was ignored in the models presented here. 



Primary sources of uncertainty in predicting fluorine nucle- 
osynthesis in AGB stars relate to the adopted reaction rates, 
especially 14 C(«,7) ls O and 19 F(a,p) 22 Ne, and the treat- 
ment of the nucleosynthesis occurring durin g the convec- 
tive t hermal pulses and the interpulse periods llLugaro et alJ 
l2004t see also Section 2.2). 



The origin of fluorine in the Milky Way 7 



Table Al. AGB 19 F and i6.17.iSq yields (M Q ). 



z 

0.0001 




X il9 F 

2.31800o-09 


X il6 

4.44786e-05 


X il7 

1.93800e-08 


1.08000c-07 


m* 
1.0 


^Tl* rem 

0.65 


Y net 19 F AGB 

0.51379E-08 


Y net 16 O AGB 

0.52869E-05 


Y net lr O AGB 

0.38546E-07 


Y net ls O A GB 

-0.59145E-08 


1.25 


0.65 


0.21202E-06 


0.69519E-04 


0.72848E-07 


-0.14220E-07 


1.75 


0.67 


0.54275E-05 


0.31623E-03 


0.60806E-06 


0.21305E-07 


2.0 


0.70 


0.10527E-04 


0.41221E-03 


0.78911E-06 


0.66983E-07 


2.25 


0.72 


0.13578E-04 


0.51952E-03 


0.41433E-06 


0.99055E-07 


2.5 


0.737 


0.45574E-05 


0.13068E-02 


0.74846E-06 


-0.17541E-07 


3.0 


0.820 


0.62016E-07 


0.13776E-02 


0.20197E-05 


-0.23173E-06 


4.0 


0.868 


0.27367E-08 


0.16630E-02 


0.36273E-05 


-0.33571E-06 


5.0 


0.924 


-0.34142E-08 


0.76699E-03 


0.19749E-05 


-0.41578E-06 


Z 

0.004 




X '19 F 

1.69218E-07 


X il6 

1.28324E-03 


X il7 

1.41477E-06 


X il8 

7.88419E-06 


m* 
1.0 


rem 

0.63 


Y net 19 F AGB 

0.94462E-09 


Y net 16 O agb 

-0.29158E-05 


Y net 17 O agb 

0.22240E-06 


Y net ls O AGB 

-0.25334E-06 


1.25 


0.64 


0.62277E-08 


-0.20669E-05 


0.49241E-06 


-0.89206E-06 


1.5 


0.646 


0.22912E-07 


0.19546E-04 


0.15093E-05 


-0.15635E-05 


1.75 


0.65 


0.17284E-06 


0.75738E-04 


0.44914E-05 


-0.25162E-05 


1.9 


0.65 


0.49608E-06 


0.95973E-04 


0.75796E-05 


-0.32065E-05 


2.25 


0.66 


0.38720E-05 


0.14749E-03 


0.11186E-04 


-0.43546E-05 


2.5 


0.68 


0.80892E-05 


0.93431E-04 


0.98989E-05 


-0.54950E-05 


3.0 


0.73 


0.68954E-05 


0.16497E-04 


0.10204E-04 


-0.69500E-05 


3.5 


0.82 


0.76727E-06 


0.62994E-04 


0.43600E-05 


-0.19736E-04 


4.0 


0.86 


0.55804E-07 


0.68820E-05 


0.17349E-05 


-0.24560E-04 


5.0 


0.91 


-0.65262E-06 


-0.19071E-02 


0.20350E-05 


-0.32029E-04 


6.0 


0.97 


-0.84234E-06 


-0.39748E-02 


0.10162E-05 


-0.39405E-04 


Z 

0.008 




Xi 19 F 

3 25423E-07 


X il6 

2 64027E-03 


X il7 

1 72075E-06 


x ns 

1 51621E-05 


m* 
1.0 


Tfl* rem 

0.5998 


Y net 19 F AGB 

0.23717E-08 


Y net 16 O agb 

-0.26083E-05 


Y net 17 O agb 

0.28610E-06 


Y net ls O AGB 

-0.39665E-06 


1.25 


0.61 


0.11401E-07 


-0.40814E-05 


0.65704E-06 


-0.15151E-05 


1.5 


0.63 


0.20168E-07 


-0.45432E-05 


0.19727E-05 


-0.28219E-05 


1.75 


0.64 


0.90125E-07 


0.22020E-04 


0.60267E-05 


-0.45889E-05 


1.9 


0.64 


0.18749E-06 


0.38094E-04 


0.95185E-05 


-0.53377E-05 


2.25 


0.65 


0.14941E-05 


0.85696E-05 


0.16553E-04 


-0.74497E-05 


2.5 


0.67 


0.30357E-05 


-0.12452E-03 


0.14976E-04 


-0.90888E-05 


3.5 


0.77 


0.22972E-05 


-0.37752E-03 


0.16132E-04 


-0.14133E-04 


4.0 


0.84 


0.66878E-06 


-0.30252E-03 


0.14747E-04 


-0.23226E-04 


5.0 


0.89 


-0.12324E-05 


-0.24193E-02 


0.68396E-05 


-0.61869E-04 


6.0 


0.95 


-0.16173E-05 


-0.57759E-02 


0.80141E-05 


-0.75982E-04 


Z 

0.02 




Xi 19 F 

4.63728E-07 


X il6 

9.60266E-03 


X il7 

3.87707E-06 


x ns a 

2.16506E-05 


m* 
1.0 


rem 

0.57309 


y net 19 F AGB 

0.36025E-08 


y net 16 O agb 

-0.61933E-06 


Y net lr O AGB 

0.34219E-06 


Y net ls O AGB 

-0.41416E-06 


1.25 


0.578 


0.15858E-07 


-0.22254E-05 


0.97603E-06 


-0.18066E-05 


1.5 


0.60 


0.25014E-07 


-0.32904E-05 


0.31910E-05 


-0.33975E-05 


1.75 


0.63 


0.29957E-07 


-0.54827E-05 


0.90611E-05 


-0.52622E-05 


1.9 


0.636 


0.28199E-07 


-0.10482E-04 


0.14418E-04 


-0.60505E-05 


2.0 


0.64 


0.27049E-07 


-0.17893E-04 


0.19332E-04 


-0.66536E-05 


2.25 


0.65 


0.12064E-06 


-0.20851E-03 


0.44894E-04 


-0.81507E-05 


3.0 


0.682 


0.34981E-05 


-0.11882E-02 


0.51877E-04 


-0.13306E-04 


3.5 


0.716 


0.49320E-05 


-0.17086E-02 


0.57427E-04 


-0.17219E-04 


4.0 


0.791 


0.14844E-05 


-0.16979E-02 


0.55128E-04 


-0.19531E-04 


5.0 


0.878 


0.86376E-06 


-0.39096E-02 


0.46522E-04 


-0.70118E-04 


6.0 


0.929 


-0.22147E-05 


-0.83059E-02 


0.78725E-04 


-0.10882E-03 


6.5 


0.964 


-0.24583E-05 


-0.97865E-02 


0.98794E-04 


-0.11879E-03 



